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Abstract 

We establish a thermodynamic Umit and Gaussian fluctuations for the 
height and surface width of the random interface formed by the deposition 
of particles on surfaces. The results hold for the standard ballistic deposi- 
tion model as well as the surface relaxation model in the off-lattice setting. 
The results are proved with the aid of general limit theorems for stabilizing 
functionals of marked Poisson point processes. 

1 Introduction 

Aggregation and growth processes associated with the deposition of particles on sur- 
faces is of considerable scientific and practical interest (Barabasi and H. E. Stanley 
[1], Cumberland and Crawford [2], Vicsek [8]). Of special interest is the morphol- 
ogy of the random interface separating the growing material from its environment. 
The surface morphology, which involves non-equilibrium growth behavior, exhibits 
scaling and is thought to evolve to a steady state without a characteristic time or 
length scale (Family [4], Family and Vicsek [5]). 

This note rigorously develops the mathematical limit theory of the random parti- 
cle deposition model in which random size l)-dimensional balls ('particles') rain 
down sequentially at random onto a d dimensional substrate Qn of volume n'^. The 
substrate Qn is assumed initially flat but it need not be oriented perpendicularly to 

Department of Mathematical Sciences, University of Durham, South Road, Durham DHl 3LE, 
England: mathew.penrose@durham.ac.uk 

Department of Mathematics, Lehigh University, Bethlehem PA 18015, USA: 
j oseph.yukich@lehigh.edu 

1 Research supported in part by NSA grant MDA904-0 1-1-0029 



1 



the particle trajectories. Assume that particle sizes are independent and identically 
distributed with diameter bounded by a constant Di. 

When a particle arrives on the existing agglomeration of deposited balls, it may 
stick to the first particle it contacts, which may result in lateral growth and 'over- 
hangs' (the ballistic deposition (BD) model). Alternatively, the particle may slip and 
roll over existing particles, undergoing vertical and horizontal displacements, stop- 
ping when it reaches a position of lower height (the surface relaxation (SR) model). 
In this model, the particle diffuses along the interface up to a finite distance, say D2 
(c.f. [4, 1]). Also, only those existing particles within distance Di + D2 of the new 
particle affect its motion. The finite displacement assumption is consistent with the 
presence of moderately sticky particles. 

Both the BD and SR models have been intensely studied in the lattice and off- 
lattice (continuum) settings (see [1] for an overview of the enormous literature). We 
focus here on the continuum setting. 

Let A' C R*^ X [0, r] denote a finite point set of particles, which are released from 
'infinity' onto the substrate R*^. Represent points in X by {X,t), where X e R'' 
denotes the spatial location (center) of the incoming particle and t its time of arrival. 
The parameter r is the 'mean deposition intensity', i.e., the average number of 
particles arriving per unit volume in R*^. 

The 'active zone' or the 'interface' of the resulting agglomeration formed through 
either the BD or SR process is defined as follows. A particle is termed 'exposed' if 
there is a positive chance that a subsequent incoming particle could strike it directly. 
The active zone is the set A{X) of exposed particles belonging to the agglomeration 
induced by X. 

Given X := {Xi,...,Xk) C R'' x [0,t], let X'-, 1 < i < k, denote the position 
of Xi after displacement. Let h{Xl) :— h{X, XI) denote the height of X^ above the 
substrate R'^. Since the maximum diameter is Di and the maximum displacement 
is D2, the height h{X) of a particle (X, t) is a random variable determined by points 
in the existing configuration within a distance + D2 of X. 

Define the height functional 




X [0,t] 



(1.1) 




(1.2) 



The squared width oi Xi e X about the mean is 



w{Xl) :=w{X,Xi) := {h{Xl) - H{X))\ XcR'^x [0,t]. 



(1.3) 
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The width functional is 



W{X):= J2 '^cM'^x[0,r] (1.4) 

and the surface width functional is 

W := W{X) := (^^) XcR'x [0, r], (1.5) 

i.e., W is the root mean square of fluctuations of heights of particles in the active 
zone. 

We are interested in the mean height and surface width functionals when X is 
the realization of a rate one homogeneous space time Poisson point process Vn,T on 
Qn X [0,t], Qn C M°'. In this case we write A{n,r) := A(Vn^r) for the active zone, 
and do not distinguish between the active zones for the BD and SR models. Let 
{Pj} denote the points in Vn,T- 

Define the random height functional H{n, r) := H{Vn,T) and the mean random 
height functional H{n,T) := H{Vn,T)- Similarly, define the random width functional 
VF(n, r) := W{Vn,T) and the random surface width functional W{n,T) := W{Vn,T)- 

Considerable research has focused on the qualitative properties of the random 
surface width W{n,T). Its analysis is rendered difficult by the presence of spatially 
correlated squared heights {h{P[) — HY, P/ e A{n, r). Both the random height and 
random width functionals involve triangular arrays of identically distributed but 
spatially correlated summands. 

Beginning with the seminal papers [4], [5], and Kardar, Parisi, and Zhang [6], 
it is widely held that the surface width exhibits both space and time scaling and is 
governed by the dynamic scaling relation 

F(n,r) ftin"/(r/n"/'^) (1.6) 

where the scaling function / satisfies 

f{x) fti x^, x<^ 1, (1.7) 

and f{x) C for a; » 1. Here, a is called the roughness exponent [1] and (3 is the 

growth exponent; a describes the dependence of the surface width in the long time 
limit r oo on the linear size n, whereas j3 reveals the dependence on time r. In 
the region t <^ n, we have 

W{n,T)^T'^ (1.8) 
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whereas for r ^ n, 

W{n,T)^n'^. (1.9) 

For both the BD and SD models there is a plethora of experimental and the- 
oretical results validating the phenomenological expressions (1.6) - (1.9). For an 
overview of the huge literature, see [1], [8]. 

Stochastic PDEs are central to the analytic treatment of continuum deposition 
models and the attendant phenomenological surface width relations (1.6) - (1.9). 
The classic linear differential equation of Edwards and Wilkinson [3], which treats 
the SR model, provides evidence for a = 1/2 and /? = 1/4 in dimension d = 1. The 
landmark paper [6] (KPZ) treats the BD model and accounts for lateral growth by 
adding non-linear terms. Subject to the validity of the KPZ methods, which are 
justified by plausibility arguments and which cannot be formally derived (see e.g. p. 
56 of [1]), one obtains in dimension d — 1 the 'exact' values a — 1/2 and P — 1/3. 
Higher dimensional values of a and f3 are not known. 

Despite the abundance of experimental results and despite the successes of PDEs, 
there is no formal, mathematically rigorous treatment of the stochastic properties 
of the surface width for the BD and SR models. 

In this note we address this lack of mathematical precision by providing a formal 
and rigorous probabilistic treatment of the mean height and surface width func- 
tionals, establishing a thermodynamic limit and Gaussian fluctuations for the mean 
height and surface width functionals. The mathematical underpinnings reside in 
the general limit theory of stabilizing functionals of marked Poisson point processes 
[7]. Although we do not obtain mathematical values for the roughness and growth 
exponents, the present approach provides additional information on a. Attention 
here is confined to the continuum, but the results also hold in the lattice setting. 

2 Statement of Results 

2.1 Thermodynamic limits 

Denote by Vr a rate one homogeneous Poisson point process on x [0, r]. Without 
loss of generality, for all n = 1, 2, ... we let Vn,T be the restriction of Vr to Qn x [0, r]. 

Intuitively, one expects that the large n limit of each summand in the random 
height functional H{n, r) :— ^p/g^(„ ^) ^(-P/) is related to a height functional defined 
on the 'infinite interface' created by the deposition of particles in Vr on the infinite 
substrate M'^. However, it is not straightforward to make rigorous sense of the notion 
of sequential deposition of an infinite number of particles over a substrate of infinite 
extent. 
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We adjoin to Vt a point {0} := {(0, i)}, with spatial coordinate the origin of 

and time coordinate uniformly distributed over [O.r]. The collection of points 
Vt U {0} forms a Palm distribution and we use it to define 'typical' height and 
squared width fluctuation functionals. Let 0' G M*^ x [0, r] be the random position 
of {0} after displacement with respect to the incoming particles Pn,T, ^(0') '■— 
h{Vn,T U {0},0') its random height, and A{n,T,Qi) the active zone generated by 

Vn^r U {0}. 

The following thermodynamic limits hold for both the BD and SR models. They 
involve notions of conditional expectation and conditional variance. Given a a ran- 
dom variable C, and an event B with non-zero probability, the conditional distribu- 
tion of C given B is the probability distribution P[{C £ P[B\] the mean and 
variance of this conditional distribution are denoted and Var[C|S]. 

Theorem 2.1 (a) For all r G (0, cxo), the random limits 

ho{r) := lim h{0') ■ l{o'eA(n,r,o)} and ^o(t) := lim l{o'eA(n,r,o)} 



exist a.s. 

(h) For allT e {0,oo), 



lim = r • EUr)] a.s. 

(c) For allr e (0, oo), 

lim H{V^,r) = = E[ho{r)\Ur) = 1] a.s. 

If we interpret h{0') ■ l{o'eA(n,T,o)} as the height of a 'typical' point in the inter- 
face, then the 'typical' squared width (height fluctuation) of a point in the growing 
interface A{n, r, 0) is given by 

W(O') • l{o'6A(n,r,0)} := w{Vn,T U {0}, O') • l{o'eA(n,r,0)}- (2.1) 

The next result shows that the typical squared width a.s. converges in the large 
n limit, and moreover the limit governs the fluctuations of the interface width. 

Theorem 2.2 (a) For all r G (0, oo), the limit 

Wo{t) := lim w{0') ■ l{0'6A(n,r,0)} 
n— >oo 

exists a.s. 
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(b) For all re (0,oo), 



n- 



lim W^{Vn,T) ■■= lim 



■oo n— >oo 



\A{n,T)\ 



1 



P' 



'/eA(n,T) 



= Var[/io(r)|eo(r) = l] 
= £;[wo(T)|eo(T) = 1] a.s. 



Remarks 



(i) Just as the random variable ho{T) gives meaning to the height of a 'typical' 
particle in an infinite interface defined by a Poisson point process, wo{t) may be 
interpreted as a 'typical' width. Similarly, E[^o{t)] is the probability that a typical 
particle belongs to the active zone. 

(ii) The scientific literature assumes the existence of the hmits lim^^oo H{n, r) 

2 

and lim„_>oo W (n, r). Theorems 2.1 and 2.2 puts this on sohd footing and identifies 
the limiting constants. In statistical terms. Theorems 2.1 (b) and 2.2 (c) say that 
the sample mean (variance) of the heights of the particles in the interface form 
a consistent estimator for the mean (variance) of the random height of a typical 
particle, conditional on its being in the interface at all, as the substrate becomes 
large. 

(iii) For both the BD and SR models, it is of great interest to estimate the con- 
stants in Theorems 2.1(c) and 2.2(b). Heuristically, one expects that E[ho{T)\^Q{T) = 
1] ^ C-T since the expected height should vary directly with the deposition intensity 
T. For large n, we expect that 0' G A(n, r, 0) with a probability inversely propor- 
tional to T, i.e., we expect E[^o{t)] ~ 1/r. Also, given io{T) = 1, one expects that 
the expectation of the squared width of a 'typical' particle in the smoothed interface, 
namely E[wo{t)\^o{t') = 1] is upper bounded by the expectation of the width of a 
'typical' particle in the non-smoothed interface. The latter is just the variance of 
a Poisson random variable with parameter r, giving E[wo{t)\^o{t) = 1] = 0{t), 
which is consistent with the simulations of Zabolitzky and Stauffer ([9], p. 1529). 
Whether this bound can be improved is the subject of widespread research. A 
rigorous mathematical treatment of E[wo{t)\^o{t) — 1] would be welcome. 

2.2 Gaussian fluctuations 

Both the height functional H and a modified version of surface width exhibit Gaus- 
sian fluctuations about the mean, as shown by the following distributional results. 
N{0, (j^) denotes a normal random variable with mean and variance cr^, and — 
denotes convergence in distribution. 
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Theorem 2.3 For all r e (0, oo) there exists :— cr^{H) > such that 



H{Vn,r)-EH{Vn,r) V 



N{0, al) 



and n-'^VeirH{Vn,T) ^ cr^- 




P/eA(n,r) 



The modified width functional Wi exhibits Gaussian fluctuations around its 
mean in the large n limit. 

Theorem 2.4 For all r e (0, oo) there exists :— cr^(VFi) > such that 



(i) If the Poisson point process Vn,T is replaced by the binomial point process 
{UnA, •••! Un^inr]}, whcrc Un,i, f^n,2, ••• are i.i.d. with the uniform distribution on Q„ x 
[0, r], then Theorems 2.1 and 2.2 continue to hold. Theorems 2.3 and 2.4 also remain 
valid, albeit with smaller variances <J^{H) and cr^(H^i). See [7] for details. 

(ii) Theorems 2.1-2.4 hold for particles of random shapes and types; the particles 
need not be balls [7]. The substrate Qn need not be a cube, but may be a more 
general region subject to a smoothness assumption on its boundary [7]. 

(iii) We have not been able to determine the dependence of <y1[H) and a'l{Wi) 
on T. The exact order of magnitude of crl{H) and 17^(1^1) could help determine the 
scaling exponents a and P in (1.6) - (1.8). 

(iv) Extending Theorems 2.1-2.4 to cases involving r and n going to infinity 
together appears beyond the scope of the methods here. Such cases are of interest 
and merit investigation, particularly since they might also help determine values for 
a and /?. 



and n-'^Y&iWx{Vn,T) 




Remarks 
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3 Methods 



To prove limit results for the mean height and width functionals, we write 

Hin, r) := := (3.1) 

and 

To prove a thermodynamic limit for H{Vn,T) it will suffice to show that the 
quotients \A{nT)\ ^^'^^^^^^ both tend to a limit a.s. The same is true for 
W{'Pn,T)- Notice that we may write 

\A{n,T)\:= J2 (3.3) 

where ^(A", x) = 1 or 0, depending upon whether x is exposed or not in the agglom- 
eration induced by X. 

To extract limits for the height and width functionals, as well as for the number 
of exposed points, we probe their behavior by adjoining {0} to A", where wc recall 
that {0} has as spatial coordinate the origin of and a time coordinate randomly 
distributed over [0,t]. The change in the height functional (the 'add one cost' [7]) 
caused by the insertion of {0} is given by 

Ah(A') H{X U {0}) - H{X) (3.4) 

with a similar definition for A^(A') and A;^^(A'). 

The displacement of the incoming particle {0} is clearly determined by the lo- 
cation and displacement of particles previously arriving within a distance D3 := 
2[Di + D2) of the origin of W^. Whether the particle ends up in the active zone 
A{X) is clearly a function of particles arriving within a distance of the origin. 

Also, the insertion of the particle {0} will only affect a later arriving particle at 
location x e R'^ if there is a sequence of points {(xj, ii)}"=i in X such that 

\xi — < I?3, and ti < ti+i, 1 < i < n — 1, 

and 

\Xn -X\ < D3. 
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When X is the homogeneous Poisson point process Vn,T, then such sequences 
exist with a probabihty decaying exponentially in \x\ [7]. Hence, adjoining {0} to 
Vn,T changes the heights h{Pl) of the re-positioned particles Pi, Pi G Vn,T, \Pi\ > 
A, with a probability which decays exponentially in A. (Thus and ^(P/) 

have exponentially decaying correlations.) The insertion of {0} thus changes the 
morphology of the active zone at distances larger than A from the origin with a 
probability decaying exponentially in A. 

We conclude that for all deposition intensities r G (0, oo), there is an a.s. finite 
random variable := Nh{t), with exponentially decaying tails, and an a.s. finite 
random variable A(r) := Ah{t) such that with probability one 

HiiVrn Qn) U{0}uA)-H{{VrnQN)^A)^ A(t) (3.5) 

for all finite A <Z {W^ \Qn) x [0,r]); see [7] for details. The 'add one cost' to 
the height functional is unaffected by changes in the particle configuration outside 
Qn X [0,r]. Following extant terminology, [7], we say that the height functional 
'stabilizes'. 

Moreover, for the same reasons, h also 'stabilizes' in the sense that there are a.s. 
finite random variables N and ho{T) such that 

h{{Vr n Qn) U {0} U A, 0') • l{o'€A((P.nQ^) u {o} u^)} = ho{r) (3.6) 

for all finite ^ C (M^ \ Qn) x [0, r]. This last fact yields Theorem 2.1(a). 

Next, given a subcube Q C M'^, let Ui^q, U2,q, ... be i.i.d. uniformly distributed 
random variables on Q x [0, r] and let Urn,Q '■= {Ui^q, f^2,Qi •••) Um,Q} be a binomial 
point process. Adjoin the particle {0} to Um,Q and consider the add one cost 

AHiUm,Q) := H{Um,Q U {O}) - //(W^.q). (3.7) 

Since H stabilizes and since Nh{t) has finite moments of all orders [7], by coupling 
binomial and Poisson point processes we can show that AH{Um,Q) has finite moments 
of all orders. In fact, for all r e [0, oo) and all p > 0, H satisfies the bounded 
moments condition 

sup sup E[Ah(W^,q)^] < oo, (3.8) 

Q me[T|Q|/2, 3t|Q|/2] 

where the first supremum runs over cubes Q C M'' containing the origin [7]. 
Likewise, the functional h satisfies the bounded second moments condition 

sup sup E[h{Um,Q U {0},0')'^] < oo. (3.9) 

Q me[T\Q\/2, 3t\Q\/2] 
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By combining the translation invariance of h, the uniform bound H(X) < 
Ci(card(A'))^, the bounded moments conditions (3.8, 3.9), and by invoking the sta- 
bihzation of h (3.6), one verifies the hypotheses of the general law of large numbers 
given by Theorem 3.2 of [7]. We thus deduce for all r e (0, oo), 

hml V h{P;)^r-E[ho{r)] a.s. (3.10) 

n— >oo fl"- ' 

P/eA(n,r) 

Similarly, using the stabilization of ^ and again using the general law of large 
numbers of [7], we may show for all r e (0, oo) that lim„_>oo ^^^^f''^ — r ■ £'[^o(''')]- 
This gives Theorem 2.1 (b). Combining this last limit with (3.10) and (3.1), we 
obtain Theorem 2.1(c). 

By invoking the stabilization of H, we may apply the general central limit theo- 
rem (Theorem 3.1) of [7] to deduce Theorem 2.3. According to Theorem 3.1 of [7], 
the limiting variance cr^iH) is non-zero provided that ^h{t) is non-degenerate. To 
check non-degeneracy of Ah{t), observe that there is positive probability that the 
set Vn,T puts no particles within distance of the origin, so that in this case the 
difference 

H{{Vr n Qn) U {0} U ^) - H{{Vr n Qn) U A) (3.11) 

is just the height of a single particle resting on the substrate. On the other hand, 
there is also positive probability that the set Vn,T is such that the addition of the 
origin results in a net increase of one exposed point to the active surface, in such a 
way that the height is increased by the height of a particle which is several layers 
removed from the substrate. This produces a different value for (3.11), showing the 
non-degeneracy of Ah{t). 

We prove Theorems 2.2 and 2.4 in a similar fashion. Define the functionals 

wi{X') := w,{X,X') := {h{X') - D{r)f 

and 

x[eA{x) 

Analogously to (3.6), the functional 

wiiVn,r U {0}, 0') := (/i(O') - Dir)f 

stabilizes in the sense that for all r G (0, oo) there is a random variable N with 
exponentially decaying tails and wq{t) a.s. finite such that 

Wi((P^ n Qn) U {0} U A, 0') ■ l{0'6A(n,r,0)} = M^) 
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for all finite ^ C (R*^ \ Qn) x [0, r]. 

Moreover, since Wi satisfies the bounded moments condition (3.8) (replace Ah 
by Awi there) and since Wi satisfies the bounded second moments condition (3.9) 
(replace h by Wi there), we obtain for all r G (0, oo), 

wo(t) = lim wi{'Pn,T U {0}, 0') • l{o'eA(n,r,o)} as- (3.12) 

n—*oo 

Since H — > D{t) a.s. by Theorem 2.1(c), it follows that 

wo(r) = lim w(K,t U {0}, O') • l{o'GA(n,T,o)} a.s. 

n—^oo 

i.e., Theorem 2.2(a) holds. 

To prove Theorem 2.2(b), we write 

Wirn,r)^ Yl {h{Pl) -H{n,T))' ^ Yl h'{P:)-\A{n,T)mn,T)r. 

P/eA(n,r) P/eA(n,r) 

(3.13) 

Dividing both sides of (3.13) by A{n, r) and appealing to Theorem 2.1(c), we see 
that to establish a thermodynamic limit for W{Vn,T) it is enough to show a thermo- 

dynamic limit for — ^^^\A{n l)\ — However, by following the proof of Theorem 2.1 
with h replaced by /i^, we find that 

Y.PleA{n,r)h\Pl) _ E[hl{T)\ 

uTo. \A{n,r)\ E[Ur)V 

Using Theorem 2.1(c) on the second term in (3.13), we obtain the first equality in 
Theorem 2.2(b). 

The second equality in Theorem 2.2(b) is just the definition of conditional vari- 
ance and it remains to show the last equality. Let En := {0' G A(n, r, 0)} for 
simplicity. We have 

{ho{T)-D{T)r 

= (lim„^ooMO')li?„-^(^))' 
= lim^^oo h\0') ■ U„ - 2D{t) lim„^oo /^(0')U„ + D{Tf 
= lim„^oo(/i(0') - D(r))2 . 1^„ + D{Tf lim„^oo Ie^ 
= lim^^oo w^i(K,r U {0}, 0') ■ 1e„ + Z^(r)2 lim^^o^ ^e^ 

by (3.12) and by Theorem 2.1(a). Thus, by the above identity 

Var[/io(r)|eo(r) = 1] = E[{ho{T) - E[ho{T)\Ur) = mUr) = 1] 
= Eiihoir) - D{T)r\ eo(r) = 1] = E[wo{t)Ut) = 1], 
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showing the last equahty in Theorem 2.2(b). 

Finally, we address the proof of Theorem 2.4. Exactly as for the height functional 
H, we can show that Wi stabilizes in the sense that (3.5) holds with H replaced by 
Wi. Invoking the stabilization of Wi and applying Theorem 3.1 of [7], we obtain 
Theorem 2.4. Non-degeneracy of Ai4/i(t) follows as in the proof of non-degeneracy 
of Aj,(t). 
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